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Abstract 

Exact bound state solutions of the Schrodinger equation for generahzed noncentral 
potential are examined by means of the Nikiforov-Uvarov method. The wavefunctions 
and the corresponding energy eigenvalues of the system are obtained analytically. The 
results examined for the potential are compatible with those obtained by the other 
methods. 
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1 Introduction 



Exact analytic studies on the Schrodinger equation for noncentral potentials have been 
considerable interests in recent years. Some of the numerous studies concerning in the 
context are available in the literature [1-7]. Two different classes of the potentials are 
considered in searching both relativistic and nonrelativistic quantum mechanical systems. 
The first case is Coulombic type potential which is used in analyzing the bound states of 
an electron within the Hartmann's ring-shaped plus Aharonov-Bohm (AB) field [8] or a 
magnetic Dirac monopole [9]. They can be used in quantum chemistry for describing the 
ring-shaped organic molecules such as cyclic polyenes and benzene as well as in nuclear 
physics for investigating the interaction between deformed pair of nuclei and the spin- 
orbit couphng for the motion of the particle in the potential fields. The second one is the 
generalized harmonic oscillator type potential known also as the ring-shaped oscillator plus 
AB-systems [10]. These are very useful for calculating the vibrational quantum levels and 
determining the hidden dynamical symmetries of quantum systems which are described 
a particle moving in the potential fields. In this study, we introduce a ring-shaped type 
generalized noncentral potential and investigate its solutions with respect to the Schrodinger 
equation. 

By replacing the Coulomb term of the Hartmann's potential with the harmonic oscillator 
one and adding to the inverse square term (or saying that it is either attractive or repulsive 
for negative and positive parameters) respectively. It also includes contribution of the 
potential term as (rcosO)"^. Thus, we define the potential in spherical coordinates as 
follows: 

V,{r, 9) = (at' + 4) + 4 (-^ + Az) > (1) 

where q is the deformation parameter [11]. In equation (1), only ^—dependent part V{d) = 
+ -^^) is often caUed the Poschl-Teller I type potential with C = ^x(x — 1), 
D — ■yA(A — 1) and Vq — [12]. Various choices of energy parameters allow us to 
investigate for some particular cases: 

For 5 = 0, the potential (1) turns into the double ring-shaped oscillator [13]; for 5 = D = 0, 
the ring-shaped oscillator [11, 13]. Moreover, for only = 0, it refers to the nonspherical 
ring-shaped oscillator [2]; for only q — {)oyC — D — the harmonic oscillator plus inverse 
square potentials, respectively. Besides these, ior B — C — D — the usual spherical 
harmonic oscillator potential is possible. 

Noncentral potentials are not usually discussed in quantum mechanical problems. The 
main reason is that most of them are not analytically solvable. In spite of this fact, some 
classes of noncentral potentials in three dimensions are solvable so long as the Schrodinger, 
Klein-Gordon and Dirac equations [14-26] with these potentials satisfy the separation of 
variables. These are also applicable to the scattering and condensed matter processes 
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[27-30]. There have been many studies involving the potentials by using the well-known 
techniques, i.e., group theoretical manner [11, 31-33], supersymmetric formalism [34-41], 
path integral method [42-48] and other algebraic approaches [49-58]. 

Throughout the study, we have focused on several purposes. One of them is to solve the 
separated Schrodinger equations for the potential (1) by using the Nikiforov-Uvarov (NU) 
method [62] analytically. Our another interest is to extract some remarkable results. 

The outhne of the paper is as follows: In Sec. II, the method is introduced briefly. In Sec. 
Ill, the wavefunction and the corresponding energy eigenvalues are obtained analytically. 
In Sec. IV, our results are reviewed and discussed. 



2 The Method: Description of the Procedure 

Let us first introduce the Nikiforov-Uvarov method. Some quantum mechanical problems 
can be solved analytically after applying the separation of variables procedure [7, 49] . Then, 
a proper transformation is chosen ) for each separated equation. Therefore, each 

of them is reduced to a generalized equation of hypergeometric type as [62] 



„, s riz) ., , aiz) , . 



a{z) 



a^z) 



(2) 



Here, f{z) is a polynomial function with the highest degree 1 and a{z) and a{z) are poly- 
nomial functions with the highest degree 2. It can also be noted that the energy appear 
as a parameter in these coefficients. Their solutions, therefore, lead to several classes of 
special functions such the classical orthogonal polynomials called as Hermite, Jacobi, Bessel 
etc. If we apply the transformation u{z) = 4>{z)y{z) to original equation (3) and arrange 
it, we have a simpler form as 



Here, one concludes 



and 



a{z)y" + r{z)y' + Xy^O. (3) 



-[\n(P{z)], (4) 



a{z) dz 



t{z) = f{z) +27i{z). (5) 

Moreover, the nontrivial solutions of equation (4) must satisfy the energy eigenvalue equa- 
tion 

An + nr' + ^[n(n - l)a"] =0, (n = 0, 1, 2 . . .) (6) 
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where r' and a" denote the first and second derivatives of them with respect to z. To 
determine 77(2;) and A = A„ by assuming 

k^\-'K\z), (7) 

and by solving the quadratic equation for 'k{z\ one gets 



' a' — t\ , / /a' — ^ 



(^^) +ka-a, (8) 



with the prime factors of a denoted as the differential at first degree. The polynomial 
expression ti{z) has to be the square of polynomials under the square root sign. The 
existence is possible only if its discriminant is zero. Thus, it takes the form without square 
root. We also point out that the determination of k is the essential point in the calculation 
of 7r{z). Our problem with a nontrivial solutions corresponds to the eigenfunctions yn{z) 
called as the Rodriguez formula 

where C„ and p{z) are the normalization constants and the weight function, respectively. 
This function must satisfy the condition 

p{z)t{z) ^ ^Jaiz)p{z)]. (10) 

Therefore, wavefunctions corresponding to discrete eigenvalues are orthonormal that they 
comply with 

j yn{z)ym{z)p{z)dz = 5mn- (H) 

If we further specify that the wavefunctions can be normalized when m — n. 



3 Calculations 

For a particle in any potential field, the wave functions ^n(r) describe the bound states and 
the corresponding energy levels E^, then the time-independent Schrodinger wave equation 
becomes 

V^^(r) + -^[E- V{r)] *(r) = 0. (12) 
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Let us first consider this equation. By putting the potential V{r)—Vq{r, 9) into the equation 
in spherical coordinates, one has 



dr 



+ 



+ 



2m 



1 d 
sin 9 89 



E 



sin 9 



89 



B 



+ 



1 8^^E 

sin^ 9 8ip^ 



1 



C 



+ 



D 



Vsin^^ cos^^ 



) |*E = 0. (13) 



By taking the total wave function as 



*(r) = ^E{r,9,(p) 



= -R{r)e{9Mcp), 
and substituting it into equation (13), one brings to the forms 



+ 



e - [Ar^ + 



2 , 7 



and 



d^Q , de 

— - + cot 61 h 

d9^ d9 



A 



C 



R{r) = 0, 



D 



sin^ 9 cos^ 9 



9 = 0, 



(14) 



(15) 



(16) 



with the well-known (/?-angle (or azimuthal) solution 

1 



m = 0,±l,±2... 



(17) 



Equations (15) and (16) are called radial and ©-angle equations, respectively and they 
are in separated forms. In equations (15) and (16), £ = e(A, S, 7) 



2mE 



A 



2mA 



B 



2mB 



, C = ^and D = and also C = C + m,j = B + A are defined. Here, 

rfi^ , andA = A(C, D, ffi) are known as separation constants. By using the NU-method, we 
seek for the solutions of the equations. 



3.1 Energy Eigenvalues 
3.1.1 Radial equation 

Let us now consider equation (15) only if £ < 0. When we apply a transformation z — r"^ 
to it, then we get the transformed hypergeometric equation 
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-Az^ - £2 - 7 R{z) = 0. 



(18) 



By comparing it with equation (2), one has 

f{z)^l, a(z)^2z, a{z)^\-Az^ -ez-i\. 
When we substitute these polynomials into equation (8), we obtain 

7r(z) = ^ ± -^4Az2 + 4(£ + 2k)z + (1 + 47). 
Prom this equation, k can be determined and we rewrite it as 



(19) 



(20) 



1 1 f (2VI^ - ^/m^), A; = -I - l^A(l + 47) 

7r(2;) = - ± - < 

2 2 [ (2VI^ + ^r+4^), A; = -1 + 1^/1(1 + 4^. 



(21) 



So the proper are the values of 'niz) that they must satisfy the condition t'{z) < and we 
have two cases: 



{{) TT = TTiiz) = I - l{2y/Az - for A; = /ci = -| - i^i(l + 47). 



In this time we get 



A = A„ = -|-^V^(l + 47)-VI 



= 2nV]4. 

As a result, it gives us the energy eigenvalue equation 



En^A 



(2n+l) + ^^l + 47 



n = 0, 1, 2.... 



(22) 



(23) 



where A = —VA. 

m 

[ii) TT = = I - 1{2VAz + VI + 47) for A; = fcs 
By recalculating A = A„ and arranging it, one gets 



+ W^(l + 47) 



En^A 



(2n + l)-l^l + 47 



n = 0, 1, 2, 



(24) 



We note that these two eigenvalue solutions are obtained by following the conditions in 
equation (21). The first one refers to the solutions of 3— dimensional harmonic oscillator, 
by taking care of the parameters given below of equation (17) as i? = and 7 = A = £{£+!). 
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3.1.2 ©-Angle Equation 

In order to solve equation (16) by the NU- method, we need to recast it into a new solvable 
form. Thus we introduce 



y"{x) + f{x)y'{x) + g{x)y{x) = 0, 



(25) 



with y{x) = v{x)p{x). Then, it is transformed to 



v"{x) + (2^ + / I v'ix) + -ip" + fp' + gp)v{x) = 0, 
VP/ P 



(26) 



where we have p{x) = exp — | / f\x)dx . By applying the above procedure for equation 
(16) and taking v ^ 6, one yields 



e"(^) + [r - (k cosec^e + d sec^e)] e{e) = o. 



(27) 



where T = (L + 1/2)^ = (A + 1/4), L is used instead of £, and k = (C - 1/4). By 
introducing a new variable as \/t — sin^, equation (27) now reads 



U-t) 1 
^"^^^ + ^0^®'^^^ ^ [t(i-t)p " " ^^^^ ^ °' ^^^^ 



Here, we have V = T/2, k = k/A and D = D/4. Again comparing it with equation (2), one 
has 



(29) 



with ( — (r + K — D^. By substituting these polynomials into equation (8), one gets 



7r{t) = -(1 - 2t) ± -^Jait^ + bit + ci. 



(30) 



Here, ai, 6i and ci parameters are equal to (4 + 16r — IGn), (4 — 16^ + 16k) and (1 + 16«;), 
respectively. Prom this equation, k is determined and we rewrite it 
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(l + 4fi;) + ^(1 + 4^) 



(31) 



^{1 + Ak)-^{1 + AD) t-^{l + AK), k=\ + \[T-{K + D)\+y{l + Ak){1 + AD). 



Proper values of Ti{t) can be chosen so that they must satisfy r'it) < 0. Two possible cases 
are valid as follows: 

(i) Having TT = 7r,(t) = i(l - 2i) - i | [^(l + 4«) + ^J(1 + AD) ] t - ^(1 + 4«;)| for 



r-{K + D)]- y{l + AK){l + AD), we get 



3 1 
-8 + 4 



1 



T-{K + D)\--^J{1 + AK)(1 + AD) 



7(1 + 4^+7(1+15) 



(32) 



By putting the potential parameters cautiously into the equation and rearranging its right 
hand sides only, one obtains 



A = 



L{L+1) 
1 



(33) 



with the compact parameter 



si = W(2ri + 1) 



(2n + l) + (y^T+4^+y(lTlD)) +- 1 + -^(1 + 4«)(1 + 4L>) 



+ {k + D). 



{ii) The other solution is 



7r = 7r2(t) = ^(l-2t)-i 



for k ^ = I + J 



{1 + Ak)-^{1 + AD) 



(1 + 4fi;)(l + AD), then we have 
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A = 



(34) 



with the compact parameter 



S2 = J(2n + l)2-(2n+l) 



1 

+ 2 



1 + Ak){1 + AD) 



+ {n + D). 



Within the framework of NU-method, the procedure which solutions are considered phys- 
ically acceptable is important. As the last two equations satisfy the rule of r'(t) < [62], 
both are acceptable solutions. 

3.2 Wavefunctions 
3.2.1 Radial wavef unction 

To determine the wavefunctions, we first consider equation (4) and we get 



(35) 



where S = (1 + a/1 + 47) and ^/z = r. From equation (10) by calculating p(z) and inserting 
it into equation (9), it stands for the generalized Laguerre polynomials as 



yn{z) = Cn Ll{z), 



(36) 



with (1 = (5 — l)/2 and the weight function is p{z) — z^^ ^^e Hence, the radial wave 
functions become 



R{z) = (f){z)yn{z) 

= z'/' e-'^^^ L^iz), (37) 

where Cn are normalization constants. When we insert the weight function p{z) and 
equation (36) in equation (11), we obtain 

C„ = / , n = 0,1,2,... 
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3.2.2 ©-Angle Wavefunction 

Prom equation (4) , one finds 



(38) 



1 + 4«) + yJ{l + AD) 



with A = 

it into equation (9), it stands for the Jacobi polynomials as 



By obtaining p{t) via equation (10) and substituting 



Z/.(t)=C',P^''^'(t), 



(39) 



where z/i = y (1 + 4/s:), t'2 = y (1 + 4D) and are normalization constants. Also, pit) 
[t{l — t)]^/'^ is used. Thus we can write 



(40) 



with = sin ^. By means of the similar procedure as in radial part, we get 



\Cn\ 



(2n + A)! 



(2n + A + 1) 



(n+f)!\ (rz!)(A!) ' 



n = 0,l,2,... 



3.2.3 Total Wavefunctions 



By following equation (14) and gathering equations (17), (37), (40) respectively, one can 
construct the total wave function in compact form 

^{z,t,^)=C^,^ z^'-'^l' e-h^^ L:{z) [t(l-t)](^+#)/^ P^' ''\t) e^'^^ (41) 

or it can also be written as 



*(r, e, ^) = Cn,n r^^-^^'^ e-^^^' L'^ir'') (sin^cos^)(^+t) P'^" "''(sin^ 6) e'"^^ , (42) 



with C„,^ = 
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4 Concluding Remarks 



In this paper, bound-state solutions of the Schrodingcr equation for (r, 6) dependent gener- 
ahzcd nonccntral potential have been investigated by the Nikiforov-Uvarov method. Both 
the wavefunctions and the corresponding energy spectra of the system have an exact and 
explicit forms. Some remarkable results are noted. By various choices of potential parame- 
ters, our results can also be reduced to the solutions of some quantum mechanical systems. 
Possible cases and results are as follows: 

(i) For example, equation (23) turns into the energy eigenvalues of the spherical har- 
monic oscillator as 



for cases A = B = C = D = and j = A = £{£ + 1), where 2n,, + £ = n (i.e., called 
as the principal quantum number). It is noted that the energy levels of the system are 
degenerate except for the ground state n — 0. The degeneracy of the system is increased by 
the larger values of the radial quantum number rir as well as the orbital angular momentum 
quantum number i. 

{ii) The deformation parameter 5, in equation (1), takes some values as [11]. If we take 
q — QorC — D — Qva. equation (1), we can thus get the solutions of harmonic oscillator 
plus inverse square potential. To run the procedure, we should follow the equations from 
(13) to (16) and consider the parameters in Section III as well as the equation (23), our 
result becomes 



{Hi) For g = 1, our results of the ring-shaped oscillator potential for [B — D — Q) and 
that of axially symmetric potential for S = are identical to those of the systems in [11, 
13], respectively. 

{iv) If all the energy parameters are chosen greater than zero, the eigenvalue equation 
is therefore written with equation (33) or equation (34) 



This is the most general energy eigenvalue result of the system. 

(f) When A = 0, in equation (1), the bound state equation (23) vanishes or collapses 
to zero. It implies the second type bound state (called as pseudo-bound state) for nonlocal 
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interactions in three-body or N-body systems [59-61]. For A ^ 0, the energy spectrum 
(23) may however approach to a continuum bound state (BCS) depending on the couphng 
parameters B, C, D and the azimuthal quantum number (m). The occurrence of a BCS 
is hence only a necessary but not a sufficient condition for such an unphysical collapse in 
few body systems. 

{vi) Equation (15) at r = 0; equations (16), (27) at = and 9 = 7r/2; and equation 
(28) at t = and t = 1 have singularities. It implies that their solutions are analytic for 

te[0, 1], and either ^^e[0, 7r/2] or ^^e[7r/2, tt] intervals. Also, the wavefunctions (41) and 

(42) vanish at points z = 0, t = and t = 1; r = and = 0, 7r/2, respectively. As 

a result, radial and angular components (r, 6) and transformation parameters {z, t) in 
wavefunctions should not take these values. 

(vii) By introducing the appropriate transformation, we get the solvable form of the 0— 
angle equation (16) which is transformed to the equation (27) and (28). Then, 6'— dependent 
part solutions of (27) refer to bound-state solutions of the Poschl-Teller potential I by 
putting r ^ E, K ^ ^xix — 1)) ^ A(A — 1) and by mapping 6 ^ x in [12]. 
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